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ABSTRACT 


The hydrodynamical-numerical prediction model developed by 
hansen is applied to the North Channel of the Amazon River 
for computation of tides and currents; the results are compared 
veh tidal prediction obtained by the harmonic method and to 
menu current measurements. A medium size grid of square mesh 
memes, 1500 m in length, represents the North Channel. The 
wing forces are the tides at the northern opening of the 
emne] near the river's mouth and the river discharge into 
the Channel at the southern end. The numerical results for 
tides were verified at three tidal stations, and it was ob- 
served that the tides predicted at the northern part of the 
channel agreed, in the worse case, within 12% of the tidal 
тапре, but those predicted at the.southern end were unsatis- 
ООШОТ т1у reproduced. This fact is attributed to the size of 
the grid which is too coarse to describe adequately the var- 
iable and irregular cross-sections and bottom topography at 
ШОО ошетип part ot the channel. The predicted currents were 


in acceptable agreement with the few available measurements. 
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I. INTRODUCTION 


A REVIEW 

Pressure gradients are the dominant mechanisms in rivers, 
narrow estuaries and some lakes which are opposed by friction 
MANE In large scale circulation the opposing forces are com- 
mations Of friction and the Coriolis force or the latter 
alone. 

Hydraulic models have been used for many years to reproduce 
am physical characteristics of rivers and have proven to be 
ρε, but space and operating problems make them a rela- 
mel) expensive tool of investigation. The appearance of high 
Speed computers made possible the numerical solution of the 
differential equations used to describe flow. A physical model 
can be transformed into a program which may be modified easily 
To Simulate different conditions and may be stored for future 
utilization and alteration, thus allowing solutions to be ob- 
tamed rapidly and at a reasonable cost. 

The hydrodynamical-numerical method developed by W. Hansen 
makes it possible to determine the motion in oceanic areas, in 
adjacent and marginal seas, in shelf areas as well as in estu- 
moles and rivers [15, 18]. In many cases the method has been 
applied successfully to the quantitative reproduction of tides 
and storm surges. The numerical treatment of the hydrodynamıcal 
equations requires the bathymetry, boundary geometry, and 


еш ото оруса! and density structure to be considered. 
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πως шосе has been applied to rivers such as, the Eider, 
KE ms anû tlie Elbe, in Germany [9, 17, 18]. The Eider is a 
tidal river in northern Germany with very irregular cross- 
Sections and a remarkably irregular depth distribution. Ap- 
plication of the HN-model with a very fine mesh, to approximate 
me cross-sections by a number of parallel canals of differing 
width and depth, produced good results with a maximum deviation 
of 30 cm between computed and observed tides, which is less 
than 10% of the tidal range. Later computations with improve- 
ments in the longitudinal direction from grid point to grid 
point made it possible to reduce the deviations between the 
computed and observed values to below 3.5 cm, which was about 
I xor rhe tidal] range [17]. Similar applications of the HN- 
model were made to the Ems and the Elbe, and the results were 
mE! cases of the same order of magnitude. For the Elbe 
the differences were smaller than 2 cm [9, 17]. 

In this study the hydrodynamical-numerical prediction model 
of Hansen is applied to the Amazon's North Channel and the 


results are compared to field observations. 


SS OBJECTIVE 

KI O b c C(Iverol this thesis is to predict by means of the 
Hansen hydrodynamical-numerical model tides and currents which 
exist in the North Channel of the Amazon River during an inter- 
Mediate river stage and to compare the numerical results with 
available field measurements. The computed tides are compared 
Pos edıction by harmonic method at three tide stations along 
the channel and the currents are compared to field observations 


selected locations. 
14 





A medium size square mesh grid is used to represent the 
North Channel's configuration and bottom topography. The 
numerical results will provide means of evaluating this rep- 
resentation in terms of grid size as well as of the merits 


and limitations of the HN-model. 


15 





ο DESCRIPTION OF THE NORTH CHANNEL OF THE AMAZON RIVER 


A. GENERAL 

The Amazon River is located in South America, where it begins 
its eastward flow from a chain of glacier-fed lakes in the Andes 
Mountains in central Peru. It emerges from the eastern foothills 
папе Andes, sweeps through the world's largest tropical rain 
Forest, and discharges into the Atlantic Ocean on the northern 
Coast of Brazil. En route to the Atlantic the river drains 
about 7 x 10° km* of land, receiving the flow of over 200 tributaries. 

The Amazon is the world's largest river in terms of drainage 
area and discharge. Its discharge is subjected to seasonal 
fluctuations. Two periods are generally considered: (1) that 
On April to June, the rain season or the high water stage; 
and (2) that from October to December, the season of low pre- 
cipitation or low water stage. 

Investigations and measurements made by the University of 
Bez, the Brazilian Navy, and the U.S. Geological Survey in 
1963 and 1964 produced some astonishing numbers for the flow 
of the river. The flow at Obidos (some 850 km from the mouth) 
during the high stage measured about 216,000 m’/sec, during low 
stage about 72,500 m’/sec, and during an intermediate stage, 
164,520 m’/sec. The average flow, estimated for the discharge 
into the Atlantic Ocean, was slightly over 175,000 m’/sec [2,19, 20]. 

The river is generally deep, and in more than ten locations 


depths greater than 100 m have been found. The width of the 
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channel ranges from about 2 km at its narrowest up to about 
10 km in some places, with long sections having intermediate 
values [2, 19]. 

Gibbs [12] and Diegues [3] have shown that the Amazon waters 
ouesd cover the Atlantic Ocean in a variable layer, 10 to 20 m 
thick, under which the oceanic water progresses towards the 
shore. Gibbs [12] showed that during the high stage of the 
river a vertical stratification ranges out over the continental 
Shelf, to distances of 80 to 230 km from the mouth, and during 
Fu low Stage it ranges out to 60 to 185 km. However, all 
vertical isohalines indicative of turbulent mixing are limited 
ντ mouth ot the river, and there is never penetration of 
Aia ter into the river as occurs in some types of estuaries 
pcs]. 

TS CHARACTERISTICS OF THE TIDES AND 

CURRENTS IN THE AMAZON RIVER 

The Amazon River seems to have the longest river tide-zone 
in the world [1]. Tides are observed up to the vicinity of 
Santarem harbor, some 800 km from the river's mouth, but during 
extremely low flows tidal effects have reached upstream to 
EENdoS. The tides are of semidiurnal type with a maximum range 
Gr. apout 5 m. Field measurements indicate that the tidal vari- 
ation decreases upstream. 

ος river flow combine to produce strong currents in 
mmy sites upstream. The North Channel (Fig. 1) has a single 
main entrance from the sea, but numerous small tributaries and 


ehanneis also offer secondary connections to the ocean. 


1 ψ 





wenan of these channels a tidal bore is a peculiar character- 
JC. 

Eraneo [8] characterizes the current in the North Channel 
mene Vicinity Of Macapa as “alternating-axial" (rectilinear 
oT TEevVerSsing), that is, the alternating current component 
superimposed on the main downstream river flow is maximum up- 
stream during the high water, maximum downstream during the 
low water and it is zero when the tide is at the mean water 
level. There is a relatively strong permanent current down- 


stream due to the enormous discharge of the river. 


STE AREA OF STUDY 

The Amazon discharges into the Atlantic in two branches. 

Hu resent Study is limited to a region of the upper branch 
which constitutes the North Channel, extending some 147 km from 
the mouth. This channel is the navigable waterway from the 
river's mouth to the ports of Macapa and Santana in the 
Duran Federal Territory of Amapa (Fig. 1). 

Mie Port ot santana is utilized mainly by the large vessels 
that carry manganese ore extracted from the mines of Serra do 
Navio, inland of Amapa. Navigation in the North Channel is 
Eu hat difficult due to natural conditions in the river, such 
as enormous floating trees, sometimes attached to floating 
Ec Urslands," pulled down from the adjacent margins by ero- 
won, strong currents, and a large tidal variation. The harbor 
santana 1s of the floating type to compensate for the large 
КЕ үле and. For the seasonal fluctuations of the river 


level. The port of Macapa consists mainly of a wharf where 
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small vessels and typical river boats load and unload passengers 
and light cargoes. 

The Channel is generally oriented in a SW-NE direction, has 
an average depth of 20 m and is never narrower than 5.5 km. In 
Some places, as in the vicinity of Macapa, the configuration of 
the river bed is quite variable, with irregular cross-sections 
and bottom topography. The bottom is composed mainly of mud 
aw sand, and the banks and bars along the river are also of 
these materials. 

A tidal bore does not manifest itself in the North Channel, 
pure the shallow sites near the mouth of the small tributar- 
шо Капа channels discharging into the main river bores are of 


very common occurrence. 
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А. 


III. THE HANSEN HYDRODYNAMICAL MODEL 


ASSUMPTIONS AND BASIC EQUATIONS 


The assumptions leading to a Hansen-type hydrodynamical 


model are the following: 


l. 
2. 


The fluid is homogeneous and incompressible; 


Pressure is hydrostatic, and thus the changes in pressure 
are due solely to changes in water surface elevation; 


Advection terms are neglectable; 


The fluid is in hydrostatic equilibrium in the vertical 
direction; and 


The geographical and vertical variations of the Coriolis 
force may be neglected. 


The basic equations for the single layer hydrodynamical 


model developed by Hansen follow from the application of these 


assumptions to the equations of conservation of momentum. 


Reference [10] presents the derivation of the equations. 


Horizontal momentum equations are integrated over depth 


lO Sive: 
b 
2 -fVe-g on χ (1) 
b 
A EU a 55 * KV*V - 10), ty) + Y (2) 
бс, S (HU) , 8 (HV) _ 
πι ху + um 0 (5) 


The wind stress components are represented by: 
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AW. [W? - μα (4) 


ШО) AW, [W + ү (5) 


The bottom stress components are given by: 


τὸ (χ) = ru(u? ν21 1/2 (6) 


ο τς ο νο s (7) 


The various terms in the above equations being defined as: 


K = coefficient of horizontal kinematic eddy viscosity 
A = drag coefficient 
WoW, = wind speed components 
ζ = surface elevation 
H = total depth (H = h + ç) 
g = acceleration of gravity 
f = Coriolis parameter 
х,у = space coordinates 
X,Y = components of external forces 
U,V = water velocity components 
1 [6 1 f5 

(U = 1f podz, V = ES 
r -uanocttomEiriction coetricient 
y Kt 


thesequdarons for bottom friction and wind stress are 
empirically developed. The wind stress terms as used in the 
model are assumed to be quadratic expressions in wind velocity, 


where X is the wind drag coefficient with a typical value of 
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3.2 x 107° [6, 7, 10]. The bottom stress is assumed to be 


nonlinearly dependent on U and V, and like the wind stress 

it was derived by empirical means. It was originally formu- 
lated for a shallow water application and the use of this 
Value for deep water is questionable, since it probably should 
be smaller in such an application. The value of r for this 


particular case was recommended by Laevastu (personal commun- 


3 3 


Mention) as 2.8 x 10 °` to 3.0 x 10 


ES FINITE DIFFERENCE EQUATIONS 

The differential equations are solved by a "leapfrog" or 
@emeral-difference scheme for time dependent solutions. At 
the boundaries the values of the velocity components and sur- 
face elevation are taken at actual points rather than from the 
surrounding points as shown in Fig. 2. 

Im Tiving forces are the tides at the open boundaries 
age the wind at the surface over the entire field. The tidal 
values are computed at each time step and introduced as new 
Values of z at each point on the boundary. 

The finite-difference form of the equations is given below. 
Using the equation of conservation of mass, the water elevation 


Me irst calculated: 


Tn) - z^ "(n,m) - f(Ht(n,m) U*(n,m) - Hi(n,m-1) U^ (n,n-1) 


* HL(n-1,m) V^(n-1,m) - Hy (n,m) V* (n,m)] (8) 


2 
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The horizontal and vertical velocity components are then 


determined from the respective horizontal and vertical momentum 


equatıons: 
U**^*(n,m) - [1 - (2rr/ HÎ” *" (n,m)) (OF (n,m)? + VF (n,m)2)¥2TF (n,m) ] 
K 
E (αμ) ] 
τε. 
* 21X (n,m) (9) 
vE*7 (n,m) - [1 - (2tr/ Hp "Un, m) (V^ (п,т) 2+ б (п,т) VF (n,m) ] 
* 
2tfÜ (n,m) - TÊ[e**" (n,m) - c^" (n*1,m)] 
+ 2 (mm) (10) 
where: 





Don) = oU (n,m) + 2 Bee one 


т (11) 


V* (n ,m) and Et (n,m) have expressions analogous to О (п,т), 


А | 
ut (n,m) - T [U* (n,m- 1) + πι + U (n,m) + U* (n*1,m)], 


κ * 
апа V* (n,m) is similarly analogous to UF (n,m). (12) 


ne time step is 21T. The total depth Шш ΟΥ A, at U or 


V points is computed as: 


BET 2T 1 СЕТ 


ү (ог у) (Nam) x h (n,m) T 2 [ζ (n,m) d a απ που) (15) 


ine effects of wind are computed by: 
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ὃν λος RS ο ο 1ο 
ΙΙΙ » ως. 205 69 
E a (μἲη η, . 1} (15) 
t H y x y Ἢ p óy? 


where P, is the barometric pressure. It is assumed that the 
pressure gradient is zero for small areas and normal conditions. 
The stability of the finite-difference scheme is governed 
ο πο Courant-Friedrich-Levy criterion which states that the 
maximum length of the time step is determined by the grid size, 


£., and maximum depth, Πο 
pc | 
Рн (16) 


пах 


At = 


For this particular grid, a value of 30 seconds was used for 
the one-half time step, At. 

When the computational area contains small regions of great 
depth, a "false-bottom'' can be used in these sections [7]. For 
example, areas deeper than 500 m may be assumed to be 500 m 
Ibis will often result in a considerable increase in 
mime Step or a decrease in the grid length if total computa- 
time 1S a critical factor. Experimental results have 
απ chat the errors introduced with this procedure are 
acceptable in some cases for practical applications of the 
model. 

The grid net is formed by three different sets of points 
as shown in Fig. 2, the water elevation (Z-points) and the 
two components of the velocity (U- and V-points). Each of 


these points has the same coordinate designation (N,M). 
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Im coastline must pass through the U- and V-points and never 
ouch the Z points. The program reads the values of the 
depth at U- and V-points as HTU and HTV, respectively. 

The Z-points are for water, land and coast designation, 
being read by the program as HTZ. They are treated symboli- 
cally as 1, 0 and -1 respectively. The HTU and HTV show the 
chart depth in centimeters, while the coastline is symbolized 
Ep Vand the land by 0, depending on their location in the 
Eid. At the open boundaries the points HTZ have symbolic 
ΙΕ of -2 and -3 for the first and second boundary when 
ШШ саре. Outside these boundaries, values of 0 for HTZ 
pr escribed. 

Me coefficient of horizontal kinematic eddy viscosity is 
related to the values of the U and V components. Hansen states 
that the computations are always stable for values of B = 
(1 «) between 0 and 0.5. This coefficient can also be used 
(ning factor: the higher the value, the smoother the 
Tent velocity distribution. In areas where the depth 
Seen bution 1S irregular, as in this case, accelerations and 
Piece irregularities appear in the model. These abnormalities 
can be corrected for by means of a proper smoothing of the 
bottom or of the water surface elevation; in this model this 
is accomplished by smoothing the water surface elevation, which 
wra numerical artifact of the model. 

їп finite-difference form, the coefficient of horizontal 
Minematic eddy viscosity is given by the relation: 


ΔΕΞ (1 - a) 
4At 


K = 


2:7 





For this model with a value of alpha equal to 0.983 recommended 
by Laevastu (personal communication), K becomes 0.46 x 10^" 
em-/sec.- 

High values of alpha must be used with caution, because 
these can allow undesirable oscillations in the model. Also, 
if too low a value is used to provide higher smoothing, the 
current velocity may be overly damped. 

As was pointed out in the assumptions of the Hansen model, 
the mean advective terms, ΙΙ have been neglected. This 
Sot necessarily a good assumption as these terms can be the 
ου σεν Of magnitude as the local acceleration terms, би./6ї. 
Mai derivation of the time averaged equations of motion, the 
turbulent contribution of the advective terms, as expressed by 
. 115 stresses, have been included using the lineariz- 
aw ddy diffusivity approximation. Hence, the non-linear 
DNE CES are not included but the amplitude effects of the 


advective terms are included by the eddy diffusivity approx- 


Matton. 


E RID SELECTION 

ΙΙ τα size is selected upon requirements of detail and 
accuracy and the availability of computer core memory and time. 
INI res not a general rule for selection of the grid mesh, 
mM Or open areas and round smooth bays, where the expected 
mimes Of current velocity and direction fields are smooth, a 
mm hesh is not necessary, but for areas where the topography 
is of primary importance, or in narrow channels, a fine mesh 


May bes required for reliable results. 


2 





τ, τισ ας model, a medium size grid was initially 
considered satisfactory (see discussion on page 6l below). 
ШЕ твае the size of the mesh will result in an increase of 
time and computer memory. A 24 x 65 unit array of square mesh 
cells having sides 1800 m in length was laid over the North 
ENannet of the Amazon River as it is represented in the Brazilian 
Netitical Chart No. 220 (BRASIL - RIO AMAZONAS, carta 220 "Da 
Barra Norte ao Porto de Santana," DHN - Rio de Janeiro). 

Bere сите performed at the Environmental Prediction 
Beer arsch Facility (ЕРКЕ) have proven that there is not much 
mierecase n accuracy for smaller grid size; for example, a 
Model for the San Diego Bay was run under three different 
meshes, and except in computation of advection and dispersion 
Gmmeolliutants, по significant differences were recorded between 


fine and coarse grids. 


ШРОТ OF TIDES AND RIVER DISCHARGE 

The tides in this model were inputed at the northern opening 
mîhe grid which was chosen and oriented such that the closest 
mE I station in that location lay over or near a grid-point, 
ms the M = 2 grid line. This tidal entrance is parallel 
to the y- coordinate. 

The tides at the station LIMAO DO CURUA (Fig. 3) were 
introduced with seven components, M2, S2, N2, K2, L2, K1 and 
MS4. Two other tide stations provide means of comparison for 
the tides predicted by the HN-model: ACAITUBA and MACAPA, 

56 km and 142 km upstream, respectively. Usually only four 


components are introduced in the model, but the values of the 
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Mev eneconstitucnts introduced in this particular case were 
considered quite Significant. 

The river discharge was inputed as permanent current, 
averaged from the bottom to the surface, in a cross-section 
near Macapa along the y-grid line at the southern end of the 
η. This input boundary was left open to insure continuity. 

There is no adequate measurement of the river discharge 
mathe North Channel at the present time. There are, however, 
Some estimates and calculations of the river discharge at 
Macapa and at the mouth and some current measurements made by 
the Brazilian Navy for navigation purposes. Table I shows a 
compilation of the more recent investigations and measurements 
for the Amazon's discharge. 

The flow of the Amazon is seasonal and can be correlated 
Ap recipitation in the entire Amazon region. Diegues [3] 
States that during the year the variation of the level of the 
mT 1S quite sufficient, the water level being closely re- 
Wied tO precipitation in the Amazon Basin to a large degree 
and the melting of snow in the Andes Mountains to a small 
ο τες. Such Variations are periodic and follow a yearly 
суете. 

The value used for this model is an estimated averaged 
pue for Macapa, 175,000 m?/sec, and it is most generally 
meeepted. Such value which is expected to represent the mean 
Sin termediate stage for the river was distributed over the 
Aros s=section in Macapa, as shown in Fig. 4, to produce the 
mean velocities (permanent downstream flow) to be inputed in 


nj o n iM € manner similar to that described by Dronkers [4]. 
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ES TUNING THE MODEL 

Tuning the model 15 а matter of patience. Several runs of 
the program were made with different combinations of alpha and 
er sborttom friction parameter. After each run both the water 
elevation and currents were compared to harmonic method and 
@ieeent measurements available. The model is left to run for 
some hours of real time with these combinations of parameters 
nu stability is attained. 

ШИ Стаб the model appears to be rather insensitive to 
Mea E rons in the bottom friction coefficient, but it responds 
weky tO small variations in alpha. The behavior of the tidal 
curve at the same location for two different values of alpha 
foo 980, a = 0.975) and two different values of bottom fric- 
RE rameter (r = 0.0028, r = 0.0030), 1s illustrated in Fig. 5. 
When changing alpha the bottom friction parameter was maintained 
Sememant and vice-versa. 

Due to a lack of adequate current measurements in the North 
ΙΙ ΠΠ] the model was mainly adjusted for the tides, and the 
Velocities of the currents were considered satisfactory when 
the tide agreement was reached and the currents approached the 


values and directions of the few available observations. 


Pee MODEL 

The calculated currents behave as described by Franco [8], 
that is, they are what is termed "alternating-axial" or some- 
times "rectilinear" or "reversing." The maximum downstream 
Telocities appear at low tide (LT) and they are greater at the 


ln cT than in the second. Similarly, maximum upstream 
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165 Occur at the high tide (HT). Slack water corresponds 
to the crossings of the mean level. 

At the grid points representing the tidal stations, the 
respective maximum and minimum predicted values of current were 
as follows: 

Mae Lando do Curua, grid Mox pot (17,2), 196 cm/sec 
IBS cm/sec for the first and second LT and 162 cm/sec and 
157 cm/sec for the first and second HT. The direction changed 
from 062° (downstream) to 242° (upstream). 

ТАШ лт: Acaituba, grid coordinates (11,26), 115 cm/sec and 
107 cm/sec for the first and second LT and 90 cm/sec and 87 
cm/sec for the first and second HT. The direction changed from 
032° (downstream) to 210° (upstream). 

(3) At Macapa, grid coordinates (7,63), 62 cm/sec and 60 
susc for the first and second LT and 42 cm/sec and 40 cm/sec 
for the first and second HT. The direction changed from 024° 
(upstream) to 201° (downstream). 

At these locations a phase lag of about 40 minutes between 
E instant of the high tide and low tide and the peak of the 
current is observed, but it is possible that this lag may be 
smaller, since the output of one hour interval of the program 
did not allow a precise determination of the exact instant of 
the low and high tide (Table II). 

τε 6 shows a plot of the current velocities at the 
three mentioned tide stations for the last 24 hours for a 


simulated 100 hours (real time) run of the program. 
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ACAITUBA 





| VEL DIR | VEL 
| ely see. | c cm/sec 

76 Js 242 80 >10 40 1964 ym 
77 i 062 43 a 42 220 

78 100 061 jan 027 33 206 

πο 165 062 60 031 J: 202 

80 196 Den 91 031 21 029 

81 192 067 109 032 43 026 

82 151 061 πι П | 51 02 

83 © 062 103 052 60 ОЕ 
84 29 242 65 πο, 62 024 

85 Wi 242 12 203 58 024 

86 162 e 71 DUM EE 41 024 

87 150 240 90 209 | 5 174 

88 98 242 84 210 32 194 

89 ЛБ 242 58 210 40 ш 
90 60 061 13 213 36 206 

91 EC 062 41 030 21 200 

92 174 0602 11 76 031 | 9 049 x 
93 183 062 97 031 4 5 029 

94 158 062 | 107 E 49 027 

95 | 98 062 103 032 x so з 
96 | 6 | 060 76 1052 | 60 | 026 

97 93 240 | 16 057 UR s OD 

| | 

98 IS A т 56 209 47 | 027 

99 | 157 242 | 86 209 00 | 043 
100 | 25  ! 194 


118 242 | 87 210 


j 


ШИШ 11 - Current velocity and direction predicted by the 

milel at the three tidal stations. The left-hand column 
ος the last 24 hrs of a 100-hr (real time) computer 

run. The right-hand columns indicate the approximate time 

ШОП Ом tides (LI) and high tides (HT). 
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iors Oi The Currents in the North Channel are shown in 
КИ ОО (a,b) through 19 in which the currents are scaled so 
Nul cm corresponds to 157 cm/sec and the arrows indicate 
the directions. The full lines on the graphs represent the 
Symbolic boundaries, the triangles represent the location of 
the tidal stations, and the ellipses the location of selected 
poants. (It should be noted that the plotting subroutines used 
for these plots are not included in the program of Appendix A.) 
Eb water elevation and the resultant current velocity at the 
Selected points are shown in Table III, which is a reproduction 
ος computer output for the special points. 

mM Fi ides predicted by the model are shown in Figs. 20 
through 22 for the last 24 hours of a simulated 100 hours (real 
fame) period. In these figures is included the difference 
between the prediction by the numerical model (HN) and the 
ПИЕ ЧАС слоп by harmonic method (HM). For Limao do Curua (Fig. 
20) the maximum difference was 19 cm which represents 4$ of 
the tidal range and a phase lag of less than 50 minutes is 
ebeerwed only at the first LT. For Acaituba (Fig. 21) the 
maximum difference was 41 cm which represents 11.7% of the 
tidal range but a phase lag of about 20 minutes is observed 
Between the HN and HM prediction. This can be partially 
Eeamained by the fact that the actual tide station lies in 
much shallower water than the closest grid point chosen to 
represent it and also because the distance from this grid 


point to the station is about 2 km. 
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BorzMacapa (Fig. 20) the results were not satisfactory. 
An attempt to improve the results, by smoothing the botton, 
mesulted in the curve labeled HNI in the figure, but further 


attempts were abandoned because they would not be realistic 


with the actual mesh size. 
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IV. DISCUSSION 


Before discussing the results obtained upon applying the 
hydrodynamical-numerical prediction model developed by Hansen 
o the Amazon's North Channel, it is worthwhile to recall some 
DENTES properties or characteristics. 

The model represents topographic features in terms of grid 
points. To obtain a good bathymetric representation a small 
mesh is required, which may sometimes lead to an extremely large 
array. It is therefore necessary to compromise between spatial 
resolution and the amount of computer core memory and time 
available. 

The time resolution is also dependent upon availability of 
µου time; in general, the smaller the time step the greater 
the computer time required. 

The flow is described by finite difference methods of solu- 
Sm of the equations of motion, which include terms of differ- 
ent orders and degrees, each one describing specific characteristics 
mile flow. These depend on the external forces applied, loca- 
Bien shape, bathymetry and properties inherent to the fluid. 

In describing the flow by these equations, it happens that some 
of the terms become dominant, so that others can be neglected 
EEUncluded in parameters such as К. τ, η Hansen's 
model neglects the advective non-linear terms, which appears 
not to cause problems if medium or coarse grids are used but 


seems to be significant for small grids, as mentioncd before. 
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The finite difference method of solution of the equations 
of motion makes continuous computation impossible. The calcu- 
lations are performed by time and spatial increments which can 
lead to enormous errors if the proper precautions are not taken. 
EN example is that of the eddy diffusivity coefficient, which 
e ected directly by the grid size and the time step since 
NS defined in terms of spatial and temporal steps. 

Шис action of viscosity in fluids is of a complicated nature. 
It may be assumed for the water motion in almost all tidal 
channels that a complete state of turbulence occurs, so that 
the flow is practically independent of the exact value of the 
Reynold's number. Because of the mathematical difficulties 
involved, the general equations of fluid motion, which include 
EE VISCOUS forces, are not generally solved in applied hydro- 
dynamics. However, the result obtained by completely neglecting 
the bottom friction is, only under the most favorable circum- 
senc es, a satisfactory approximation to the actual fluid 
motion. Therefore, to determine the flow conditions in a 
Sel, the influence of friction cannot be neglected [4]. 

Wo include friction, because the velocity varies in the 
οπές] as well as in the vertical directions in the cross- 
sections of the channel, there remains only the possibility 
of experimentally determining it by tuning the model so that 
predicted values match as closely as possible the observations. 

In tidal channels and inland rivers, the river beds are 
usually irregular and form hydraulically rouph Suit aces wera - 


bulence being generated by the irregularities of the bottom. 


5:9 





At the bottom the velocity is zero and increases in the 
vertical direction until it reaches a maximum value at some 
Get ance from the free surface. From this point to the sur- 
ET the velocities decrease slightly. Observations in very 
wide channels have shown that the sides of the channel ‘have 
Meter rcally no influence on the velocity distribution in the 
Bemteral region [4]. 

In engineering and practical problems the determination 
of the mean velocity in a cross-section is usually sufficient, 
thus making it necessary to have a relationship between the 
mean velocity and the frictional forces. 

Hansen's model assumes that the bottom stress is directly 
ELcendent on U and V but only indirectly on X and t, and the 
same constant friction coefficient is used in all computations 
(for this case r = 0.0029). 

Emu the setting of the boundaries and the representa- 
mot the coast line and margins is also a matter of personal 
mm nent, and because of this, often the representation of the 
real boundaries along the grid-points is somewhat difficult and 
may be questionable. 

For cur particular case, some facts have contributed to 
make adequate conclusions impossible. First, the river dis- 
charge is inputed as a permanent flow using an estimated value, 
because at the present time there are no available measurements 
of the discharge in the Amazon's North Channel. The distribu- 
tion of the mean velocities in the cross-section was made at 
only five grid-points which received relatively large values, 
il Ml, atfecting the resolution. 
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Second, there is insufficient adequate current measure- 
ments to provide a better adjustment of the predicted current 
Welocities, either for the improvement of the estimated river 
discharge, neither for posterior verifications. 

Third, the only tidal station in the middle portion of the 
Channel, ACAITUBA, had its constituents determined in a very 
short period of time, and the harmonic prediction for that 
tion, using only few constituents, certainly is not the 
best. Furthermore, as mentioned before, this station is lo- 
cated about 2 km from the closest grid-point, which in turn 
Exocated in much deeper water. 

In spite of these facts the model produced numerical results 
mete, if not optimal, are at least acceptable. 

The tides predicted at the northern station, LIMAO DO CURUA, 
agreed with the harmonic prediction within 4$ of the tidal range 
and no phase lag is observed. The maximum difference for the 
Station at ACAITUBA is about 12% of the tidal range, and a 
r lag of about 20 minutes is observed, but both may be 
attributed to the fact that the actual station is over shallower 
water than it is represented in the grid-point and to the dis- 
mimee from the grid-point to it. 

At MACAPA, however, the results are uns atistactor Es 
δρ at least partially, attributed to the fact that the 
Prd size is too coarse to describe correctly the bottom to: 
pography and the cross-sections in the southern 18km of tire 
Channel, both very irregular. In fact, smoothing the Taree 


eradients between the depths at U- and V- grid-points of same 
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coordınate designation by making them shallower resulted in a 
Mater curve, as is illustrated in Fig. 22. No further im- 
provement was tried because, due to the actual mesh size, the 
new "false-bottom" would not be realistic, and the procedure 
was abandoned. A similar situation is reported by Ramming [17] 
when modelling the Eider River in Germany: the first attempts 
mamreproduce the tides in areas of very irregular cross-sections 
Budedepths were not satisfactory, in particular with regard to 
mae steep flood part of the tidal curve. Differences of about 
БОКЕ were recorded in a 3.8 m tidal range. Representing the 
Exu Section by a fine structure of depth distribution, 1i.e., 
WME a finer grid, led to results in which the maximum devia- 
Bon was less than 10% of the tidal range and later improvements 
reduced this deviation to less than 1% of the tidal range [17]. 
MMesprediction of currents seems reasonable. At least at 
EE four points where our current predictions were compared 
wene tual observations made by the Brazilian Navy, in identical 


E lions of river stage, the results are quite satisfactory. 
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V. SUMMARY 


The hydrodynamical-numerical prediction model developed by 
ansen was applied to the North Channel of the Amazon River 
prediction of tides and currents. 

mie tides predicted by the model agreed, in the worse case, 
pomuthin 125 of the tidal variations at two stations, but the 
mus at a third location weren't satisfactory. This is 
attributed to the fact that the irregular shape of the channel 
in the vicinity of this station was not adequately represented 
by the grid. Improvement in the results was obtained by smooth- 
E the bottom at the U- and V- grid-points in order to decrease 
Mme gradients existing between these points, but no further 
adjustments were tried because they would be not realistic with 
the mesh size used. The currents predicted by the model appear 
ШО satisfactorily reproduced, but more measurements are 
meauared to allow further conclusions. 

EN ἔπε intention of the author to continue the investi- 
Patrons in this field to achieve improved results, both by 
decreasing the mesh dimensions and by using new measurements 
to be obtained in the near future. Continuation of the study 
will also be in the form of a two-layer HN-model which the 
author is preparing for the Amazon Estuary, for which more 


data and information are available. 
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